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ABSTRACT 

We present optical and infrared characterization of the polluted DAZ white 
dwarf GALEXJ193156. 8+011745. Imaging and spectroscopy from the ultravio- 
let to the thermal infrared indicates that the white dwarf hosts excess infrared 
emission consistent with the presence of an orbiting dusty debris disk. In addi- 
tion to the five elements previously identified, our optical echelle spectroscopy 
reveals chromium and manganese and enables restrictive upper limits on several 
other elements. Synthesis of all detections and upper limits suggests that the 
white dwarf has accreted a differentiated parent body. We compare the inferred 
bulk elemental composition of the accreted parent body to expectations for the 
bulk composition of an Earth-like planet stripped of its crust and mantle and find 
relatively good agreement. At least two processes could be important in shaping 
the final bulk elemental composition of rocky bodies during the late phases of 
stellar evolution: irradiation and interaction with the dense stellar wind. 

Subject headings: circumstellar matter — planet-star interactions — stars: abun- 
dances — stars: individual (GALEX J193156. 8+011745) — white dwarfs 



Introduction 



Recent studies have solidified the paradigm that metal-ri ch white dwarfs are po l luted 
by remnant rocky bodies from their planetary systems (e.g., IZuckerman et a/.l 120071 : lJura 
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_ _ _ Farihi et all boool . boiOal lbl: lOufour et all I2OI0I : iKlein et all I2OI0I : iMelis et all I2OI0I : 
Zuckerman et al. 201ol . and references therein). It is now desirable to understand the mass, 
nature, and history of the accreted parent body at each white dwarf. The larger the degree 
of pollution, the greater potential for insight into extrasolar rocky bodies. 



Venues et all (120101 ) describe a heavily polluted DAZ white dwarf - GALEX J193156.8+011745 
(hereafter GALEX1931) — exhibiting atmospheric magnesium, silicon, oxygen, calcium, and 
iron (Mg, Si, O, Ca, a nd Fe respectively). Additionally, Two Micron All Sky Survey (2MASS 
Skrutskie et allUOOw photometry suggested excess near-infrared emissi on. Abundance anal - 
ysis and the shape of the putative infrared excess was interpreted by IVennes et all (I2OIOI ) 
as consistent with a mid L-type dwarf that orbits and pollutes the white dwarf. However, 
near-infrared emission fr om a prominent dust disk is also consiste nt with the spectral energy 
distribution presented by I Venues et all (I2OIOI ). iDebes et all (I2OIOI ) report Wide- field Infrared 
Survey Explorer (WISE) imaging and optical spectroscopic monitoring of GALEX1931 that 
support a disk interpretation. 

Here we present optical through infrared imaging and spectroscopy of GALEX1931 
designed to better constrain the presence of infrared excess emission and to probe the material 
polluting the stellar atmosphere. 



2. Observations 



2.1. Nickel Optical Imaging 



Optical imaging was performed on UT 25 May 2010 at Lick Observatory with the 40" 
Nickel telescope. These observations used the facility's Direct Imaging Camera (CCD-C2), 
a 2048 X 2048 pixel detector with 15 fim pixels. The 0.184" pixel"^ plate scale affords a 
field of view of ~6.3' squared. The detector was binned by two in rows and columns and 
was readout in fast mode. Occasional clouds passed through the region of the sky being 
observed; frames affected by clouds were not used in the final analysis. 



GALEX1931 was observed in each of the BVRI filters f lBesselllll990r ). For each filter a 4- 
step dither pattern with 10" steps was repeated with 20 second integrations per step position. 
The flux calibrator star SA 111-1925 (ILandoltl 120091 ) was visited following the completion of 
a dither sequence for each filter. Images are reduced by median-combining all good frames 
to obtain a "sky" frame and subtracting this sky frame from each image. Sky-subtracted 
images are then divided by fiat field frames obtained by imaging the twilight sky (BVR) 
or the illuminated telescope dome (I). Each science frame is registered using bright stars in 
the field and then all science frames are median combined to yield the final reduced image. 
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Detector counts for GALEX1931 and SA 111-1925 are extracted with an aperture that yields 
~85% encircled energy (with a negligible correction between the two stars). This is achieved 
by extracting counts for both stars with a 6 pixel (2.2") radius circular aperture. Sky is 
sampled with an annulus extending from 20-100 pixels. Uncertainties are derived from the 
dispersion of measurements made from individually reduced frames. The uncertainties for 
SA 111-1925 are propagated into the final quoted uncertainty for GALEX1931. The final 
V-band measureme nt from the Nickel observations i s less accurate than that reported in 



Vennes et al\ ( 120101 ): we thus adopt the IVennes et al\ ( 120101 ) V-band measurement. 



Nickel images are presented in Figure [T] and p hotometry is repor ted in Table [T] which 
also includes Galaxy Evolution Explorer (GALEX; iMartin et a/.ll2005[ ) fluxes. 



2.2. ISAAC Near-Infrared Imaging 



Near-infrared imaging was performed with ISAAC (iMoorwood et a/.lll998l ) mounted on 
the VLT-UT3 telescope. Data were obtained on UT 29 October 2010 at Cerro Paranal in 
excellent conditions, resulting in stellar images with FWHM <0.5" in the Jg- and H-bands 
and <0.4" at Kg. The Jg filter was employed due to a red (K-band) leak in the standard 
J filter. Utilizing individual exposure times of 10 s, a 9-point dither pattern with non- 
uniformly spaced posit ions inside a 30" boun ding box was executed at each of JsHKg. The 
standard star EG 141 (IHawarden et a/.ll200ll ) was observed in an identical manner but with 
5 s individual exposures. 

These data were reduced at the telescope by the Gasgano ISAAC pipeline. A sky frame 
is generated from the median of all the science images and subtracted from individual ex- 
posures. A master fiat field is employed to correct pixel to pixel gain variation. Individual 
science frames are registered using 2D cross- correlation and averaged to yield a fully reduced 
image. Photometry is performed on both GALEX1931 and the standard star with the IRAF 
task apphot using an 8 pixel (1.2") aperture radius and 20-30 pixel sky annulus. Flux mea- 
surements for GALEX1931 are corrected for a mild airmass difference between GALEX1931 
and EG 141 that amounts to <0.02 mag, and fiux-calibrated on the Vega system. Although 
S/N>300 was achieved for the science target at all three bandpasses, a conservative error 
of 5% is chosen for the derived magnitudes to compensate for uncharacterized filter trans- 
formation furiction s for the ISAAC filters and other uncharacterized systematics (see e.g., 
Leggett et aZIbooel ). 



Imaging photometry of GALEX1931 was also performed at L' using both telescope 
nodding and secondary mirror chopping. An ABBA nodding pattern was repeated with 
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several arcseconds between each nod position and approximately 30 s spent at each nod step 
(or 120 s for the entire cycle). At each nod position, chop throws of 20" were utilized at a 
frequency of 0.43 Hz with individual exposure times of 0.11 s and 9 coadds at each mirror 
position. Additionally, the position of the ABBA pattern on the array was non-uniformly 
dithered between cycles within a 10" box. The total on-source time for GALEX1931 was 1188 
s, consisting of 30 choppi ng cycles nested wit hin 10 nodding cycles of 4 positions each. The 



L' standard star FS 148 (ILeggett et a/.l 120031 ) was observed similarly with a total on-source 



time of 118.8 s over 1 complete nodding cycle. 

A combined and reduced L' frame was generated at the telescope by the pipeline. The 
two coadded and chop-subtracted images produced at each consecutive pair of nod positions 
are combined, and then divided by a master flat fleld. Image alignment is achieved by 
compensating for the executed telescope offsets between nodding, and the two images are 
averaged. The difference in chopped images effectively removes the sky background, and the 
observing pattern results in one positive and two negative stellar images in the recombined 
frame. A flnal, reduced image is constructed by combining (aligning and averaging) all the 
frames obtained in this manner. Photometry on both GALEX1931 and the standard star is 
performed with apphot using an aperture radius of 20 pixels (1.4") and sky annulus of 30-50 
pixels. The difference in airmass between GALEX1931 and FS 148 implies an extinction 
correction of less than 0.01 mag, and the flux of GALEX1931 is calibrated on the Vega 
system. A total photometric error is calculated by the quadrature sum of an (assumed) 5% 
calibration uncertainty and the measured S/N = 9.2 at the photometric aperture radius. 

ISAAC photometric results are listed in Table [H Figure [1] shows the ISAAC images of 
GALEX1931. It is noted that the L'-band flux detected from this imaging differs signiflcantly 
from the 3.35 fl ux detected toward GALEX1931 by the WISE satellite of 1.02±0.04 mJy 



(jPebes et a/.ll2010l ). Possible reasons for this discrepancy are discussed Section [321 



2.3. Optical Spectroscopy 



Grating resolution optical spectroscopy was performed at Lick Observatory with the 
KAST Double Spectrograph mounted on the Shane 3-m telescope. Table [2] lists the obser- 
vation date and instrumental setup. Data are reduced using standard IRAF longslit tasks. 
Absolute flux calibration of th e KAST data is acco mplished by scaling to the Nickel B- and 
R-band measurements and the lVennes et al\ (120101 ) V-band measurement. 



Keck HIRES fIVogt et a/.l 119941 ) optical echelle spectra were obtained for GALEX1931. 
Table [2] lists the observation date and HIRES instrumental setup. Data were reduced using 
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the MAKEE software package which outputs hehocentric velocity corrected spectra shifted 
to vacuum wavelengths. After reduction and extraction, high order polynomials are fit to 
each order to bring overlapping order segments into agreement before combining both HIRES 
exposures and all orders. 



2.4. FIRE Near-Infrared Spectroscopy 



Near-infrared spectro scopy was obtained with the Folded-port Infrared Echellete (FIRE; 
Simcoe et al\ |2008| . |2010| ) mounted on the Magellan 6.5-m Baade telescope. Observation 
parameters are listed in Table [2J 

Prism-mode observations were performed with a 0.6" slit aligned with the parallactic 
angle. One ABBA nod-pattern was obtained for GALEX1931. A single nod pair of 1 s inte- 
gration per nod was obtained for the telluric calibrator star HD 189920 (AO V). HD 189920 
was too bright to be observed directly through the slit. The telescope was offset by 6" from 
HD 189920 and defocused until un-saturated spectra could be obtained. Although such an 
observation strategy would prevent accurate absolute flux calibration, the r elativ e flux cal- 
ibration should not be impaired. Data reduction followed iBurgasser et al\ (120101 ) with the 
exception of the dispersion solution being fit with an 8th-o rder Legendre polynom ial (essen- 
tially equivalent to the 5th-order cubic spline employed by IBurgasser et a/.ll2010l ). Absolute 
flux calibration of the GALEX1931 prism data is accomplished by scaling its spectrum to 
the ISAAC JsHKg-band measurements. The FIRE prism data are presented in Figure [31 



3. Results and Modeling 
3.1. Optical Spectroscopy 



Parameters for the white dwarf star are derived from model fits to all Balmer lines in 
the blue-side KAST data (Fi g ure E l). The method relies on the "spectroscopic technique" 
developed by [Bergeron et al. fll992h and described at length in iLiebert et al\ (120051 ) and 
references therein. We find that the KAST data are best reproduced by a white dwarf 
with parameters Teg = 23470±300 K and log^f = 7.99±0.05 where g has units of cm s^^. 
We note that GALEX1931's effecti ve temperature obta ined from our fitting method differs 
significantly from that reported in IVennes et al\ (120101 ) — Teg = 20890±120 K; the origin 
of this discrepancy is not clear at this time. By matching a model white dwarf atmosphere 
with our derived parameters to the optical spectra and photometry we derive a distance to 
the white dwarf of 56±3 pc — this agrees within the quoted uncertainties with the value 
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calculated by IVennes et al.\ ( 120101 ) . It is noted that th e spectra and photom etry also favor a 
hotter Teff for GALEX1931 t han the value re p orted in lVennes et all (|2010| ). Cooling models 
similar to those described in Fontaine et all (l200l[ ). except with carbon-oxygen cores, are 
used to estimate a cooling age of 30 Myr. The white dwarf mass and enve lope mass are 
estimated to be Mwd = 0.63 M0 and logg = —16.1 (where q= Menv/MwD; iKoesterl |2009| . 
D. Koester 2010 private communication). 

With these parameters for GALEX1931 we proceed in fitting the metallic absorption 
lines detected in the HIRES spectrum. We use a lo cal thermodynamic equilli brium (LTE) 
model atmosphere code similar to that described in iDufour et al\ ( l2005l . 120071 ). Absorption 
line data are taken from the Vienna Atomic Line Databas^. Using the effective temperature 
and surface gravity determined from fitting the Balmer lines we calculate grids of synthetic 
spectra for each element of interest. The grids cover a range of abundances typically from 
log[?7-(Z)/n(II)]= —3.0 to —7.0 in steps of 0.5 dex. We then determine the abundance of 
ea ch element by fitting the various observed lines using a similar method to that described 



m iDufour et all ( l2005l ). Briefiy, this is done by minimizing the value of which is taken 
to be the sum of the difference between the normalized observed and model fiuxes over the 
frequency range of interest with all frequency points being given an equal weight. This is 
done individually for each line and the final adopted abundances (see Table [3]) are taken to 
be the average of all the measurements made for a given element after removing outliers. 
Uncertainties are taken to be the dispersion among abundance values used in the average. 
If this dispersion is less than 0.10 dex, then the abundance uncertainty is set to 0.10 dex. 



As in lVennes et all ( I2OIOI ). we find contributions from Ca, Mg, Si, and Fe. The HIRES 
spectra do not exte nd sufficiently into the red to confirm the presence of O, so we use the 
Vennes et all ( I2OIOI ) O I A7777 triplet equivalent widths in place of our own measurements 
to estimate the O abundance. We make synthetic spectra for the O I triplet in the same 
way as for the ot her elements. The equivalent widths of the synthetic O I lines are then 
compared to the IVennes et al\ ( I2OIOI ) measured values to derive the abundance reported 
in Table [31 Abundance naeasur ements for Ca, Mg, Si, and O are in agreement with those 
reported in IVennes et al\ (120 lOf ) not withstanding the significant T^s di fference between our 
white dwarf parameters and theirs. The revised abundances quoted in IVennes et all (1201 if ) 
are als o in agreemen t with our own, except for Mg which is less abundant by ~0.3 dex due 
to the IVennes et al\ (120111 ) inclusion of pressure broadening and non-local the r modynamic 
equilibrium effects. However, we find an Fe abundance higher than I Vennes et all (j2010ll201ll ) 
by ^0.4 dex. It is noted that fits to the Fe lines assuming a T^s of 21,000 K yield an Fe 



^http : //vald . astro . univie . ac . at/~vald/php/vald .php 
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abundance consistent with that reported in IVennes et al\ (120101 ) . New elements discovered 
in the HIRES data are chromium identified from two Cr II lines and manganese identified 
from one Mn II line (Figure H]). Aluminum may potentially be present in the HIRES data as 
well (Figure H]), however the S/N is marginal so we report Al only as an upper limit (Table 
[3]). Abundances for all elements identified to date in the atmosphere of GALEX1931 are 
reported in Table [3] as well as upper limits for other elements of interest (see Section Hj). 

It is noted that the GALEX NUV and FUV measurements appear discrepant with the 
white dwarf parameters derived above. We find that t he model can be broug ht into agreement 
with the data by reddening the white dwarf with the ICardelli et al\ (Il989[ ) extinction curve 
assuming E(B— V)= 0.03 and interstellar reddening. However, such reddening places the 
model at the limit of the 3a uncertainties quoted for the optical photometry. It is likely that 
the true color excess toward GALEX1931 is less than 0.03 and that some combination of 
reddening and metal line blanketing (from elements detected and as of yet undetected) in 
the UV-spectrum of GALEX1931 are responsible for the depressed GALEX measurements 
relative to the model. 



3.2. Infrared Excess Modeling 



Table [T] lists measured ultraviolet, optical, and near-infrared fiuxes for GALEX1931. 
These data are presented with the KAST and FIRE prism- mode data in Figure [31 Compari- 
son of the Table [T]J.HK, magnitudes t o the 2MASS magnitudes (J=14.66±0.05, H=14.55±0.09, 
Ks=14.45±0.10; 



Skrutskie et al. hood ) suggests that the 2MASS measurements for GALEX1931 



suffer from contamination in the photometric apertures and/or annuli used to calculate the 
sky background and hence are inaccurate. Looking at Figure [H one can identify at least 
two sources that may have fallen within the 2MASS ap ertures in th i s croAv ded field near the 
Galactic plane. Similarly, the WISE data reported in iDebes et al\ (120101 ) could have been 
contaminated by these same sources. Due to the discrepancy between the WISE data and 
the high angular resolution data presented herein, we do not include the WISE data in our 
analysis. The ISAAC Jg- and H-band measurements agree to within their errors with the 
white dwarf photospheric model while the Kg-band measurement confirms the upward slope 
of the FIRE spectrum demonstrating ex cess emission t hat b egins at wavelengths near 2 ^m. 
This result challenges the conclusions of lVennes et al\ (120101 ) whose suggestion of an L5-type 
dwarf companion relied on the H— Kg color of the 2MASS photometry. Figure illustrates 
that the excess emission implied by the contaminated 2MASS photometry is potentially 
consistent with an L-type companion. The ISAAC photometry, which is of higher fidelity, 
clearly indicates that a roughly mid L-type object cannot be the cause of the excess infrared 
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emiss ion. Later-type companions are ruled out by the ISAAC photometry (e.g.. iFarihi et al. 
20051 ) and FIRE spectra: the spectra smoothly connect the ISAAC JgHKs data points and 
show no evidence for the e xpected spectral features arising from a sub-stellar companion 
dParihi fc Christopherl [20041 ). 



Given the volume o f previous results linking white dwarf metal pollution and dusty disks 
(e.g., iFarihi et a/.l 120091 ) . it seems reasonable to conclude that the observed excess infrared 
emission detected toward GALEX1931 emanates from a disk that orbits (and is accreted by) 
the white dwarf. A model for the photospheric emission of GALEX1931 is computed with 
the parameters described in Section 13.11 and including the detected metals with abundances 
as listed in Table [31 To fit the observed near-infrared spectra and photometry, we a dd to the 
white dwarf model an optically thick, fiat dust disk f lJural l2003l : lJura et a/.ll2007al ) with an 
inner disk temperature (Tinner) of 1400 K, outer disk temperature (Toutcr) of 1200 K, and an 
inclination angle {€) of 70° (where i = 90° is edge on to our line of sight). In this fit we set 
Rwd/D = 5.3x10^^^ (see Section |3TT]) . This disk and star model is shown in Figure [3] and 
reproduces well the observed data. It is noted that the fit presented here is not unique and 
that the temperature of the outer region of GALEX1931's disk is weakly constrained by the 
present data. 



For the case of optically thick, fiat dust disks Equation 1 from lJural (|2003[ ) can be 
used to determine at what radial distance from a white dwarf dust particles of certain 
temperatures reside. The inner disk temperature of 1400 K suggests particles residing as 
close as 25 R^. for a flat disk geometry. A temperature of 1400 K is close to the sublimation 
temperature for some silicates, the expected mineralogic al constituent of GALEX1931's dusty 
disk given the stoichiometry of other white dwarf disks (jjura et a/.ll2009al : iReach et a/.ll2009l ) 
and the metals detected in its atmosphere (Table |3]). Thus, the modeled inner disk radius 
could be physically motivated by sublimation of grains that will eventually be accreted 
by GALEX1931. Our data set is not comprehensive enough to rule out a warped disk 
morphology which could increase t he inner disk radius while still reproducing the observed 
near- infrared excess emission (e.g.. I Jura et a/.ll2007bl . l2009al . and references therein) . It is 



noted that if a fiat dust disk geometry accurately describes the distribution of GALEX1931's 
orbiting debris, then the disk is likely narrow as fiat dust disk models with Toutcr^ 1000 K 
(Router< 40 Rwd) are not consistent with the data. Su ch a result suggests GALEX1931 may 



host a narrow distribution of disk material (e.g., see IFarihi et a/.ll2010bl ). Observations at 



longer wavelengths should solidify which of the above interpretations are valid. 



In their red echelle spectra of GALEX1931. IVennes et all ( l2010l ) do not identify double 
peaked emission lines from the Ca II in frared triplet indicative of an orbiting g;aseous disk 
( iGansicke et al\ l2006l : iMelis et al\ l2010l ). Modeling and results presented by iMelis et al. 
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( 120101 ) and lJural (120081 ) suggest that gas disks around white dwarfs are fed by an influx of 
rocky objects into the Roche hmit of a white dwarf star. The absence of such emission 
hues and hence an obvious gas disk orbiting GALEX1931 are suggestive of its dust disk 
originating in the tidal shredding of one large rocky object. 



Discussion 



Modeling of the ultraviolet through infrared spectral energy distribution of GALEX1931 
provides strong evidence that an orbiting dusty circumstellar disk emits the detected excess 
infrared emission. Following earlier results, we assume that this dusty disk has its origin in 
the tidal disruption of a rocky object that breached GALEX1931's Roche limit. The ultimate 
fate of such material is to be accreted by the white dwarf star. Under this assumption we 
can deduce the bulk elemental composition of such an object from the metal abundances 
reported in Table El To properly convert the measured abundances to the abundances 
at the time of accretion we must include the effects of diffusion, radiative levitation, and 



the relevant accretion phase (IKoesterl l2009l : IChayer fc Dupuid l2010l : lJura et al\ l2009bl . and 
references therein). The presence of a disk orbiting GALEX1931 and the extremely short 
elemental diffusion constants (Table [3]) imply that we are observing elemental abundances 
during the steady-state accretion phase. In this phase, accretion of metals has reached an 
equilibrium with the rate at which th ese metals diffuse out of the white dwarf's photosphere. 
The work of IChayer fc Dupuid (120101 ) suggests that for accretion rates >10^ g s~^ the effects 
of radiative levitation are negligible — all measured accretion rates and upper limits exceed 
this threshold value (Table |3]) and hence we ignore radiative lev itation. Th u s, we d etermine 
the accreted abunda nce ratios by em ploying Equation 16 from lJura et al\ (l2009bl ) and the 
diffusion constants of iKoesten (120091 ) computed specifically for the GALEX1931 parameters 
determined in Section l3?n (D. Koester 2010, private communication). Accreted abundances 
relative to Fe are reported in Table I3l 

We desire to make a rough estimate of the mass accreted to date onto GALEX1931. 
Accretion rates for each element are computed under the assumption that material accretes in 
the steady-state; these values are reported in Table I3l The total accretion rate summed over 
all elements is ~3.9 x 10^ g s~^. An estimate of how long GALEX1931 has been accreting 
from its dusty disk comes from the percentage of DA stars with effective temperatures similar 
to GALEX1931 and atmospheric pollution. It is assumed that this percentage is linked to the 
length of time in which such DA white dwarf stars would show metal pol lutio n and hence be 



pollut ed from a surrounding disk of material. From iKoester et al\ (|2005[ ) and iKoester et al. 



( 120091 ) we find that 5 out of 116 or 4.31% of DA white dwarfs within the temperature range 
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of 19,500 K < Teff < 25,000 K are polluted with metals (including GALEX1931). Taking the 
cooling time of a white dwarf in the above T^s range to be equal to that of GALEX1931 (30 
Myr, see Section [XT]) , and assuming all white dwarfs accrete from a disk during this interval, 
it is determined that the disk hfetime for GALEX1931 is 1.3 Myr. Assuming GALEX1931 
has been accreting at a constant rate for ~1 Myr, we estimate that ~2 x 10^'^ g of material 
has accreted onto the white dwarf star. This is close to the second most massive asteroid in 
our Solar System, Vesta (~3 x 10^'^ g). It is noted that numerous assumptions are included 
in this estimate and that it could increase or decrease significantly and that the estimate 
does not include the mass remaining in the disk. 



Following iKlein et all ( l2010l ). we attempt to balance the accreted oxygen abundances 
under the expectation that all oxygen was contained in rocky oxygenated minerals (e.g., 
MgO, Al.O.s, Si O,, CaO, TiOa, CrsOg, MnO, FeO, FeaOg, and NiO) and water (H2O; see 
Jura fc Xull2010l . for a discussion of water survival in rocky ob ject s during pos t -main sequence 
evolution). The criterion that needs to be met is Equation 3 of iKlein et al\ (120101 ): 



bal 



^ QjZ) n{Z) ^ ^ 



where element Z is contained in molecule Zg(^z)Op(^z)- In this equation one must consider 
only hydrogen accreted from the parent body (Hacc)- We find that the oxygen accreting onto 
GALEX1931 can be accounted for almost entirely by inclusion in rocky minerals. There is, 
however, the ambiguity of how Fe was partitioned in the parent body that accreted onto 
GALEX1931. In our mo del Fe could have been contained in any one of metallic Fe, FeO, or 
Fe203 (see discussion in lKlein et a/.ll2010l ). Depending on how we divide Fe into these three 
states, there can be anything from a significant deficit in the oxygen balance {Obai = 0.58 
instead of 1.0, indicating too much oxygen for the number of heavier element atoms) to a 
significant overestimate of rocky minerals {Obai = 1-33 instead of 1.0, indicating tog little 
oxygen for the number of heavier element atoms). Assuming as did lKlein et al\ (I2OIOI ) that 
half of the Fe is in the form of metallic Fe and half in the form of FeO, we arrive at a value 
of O;,a;=0.83. With this value we constrain the amount of water present in the object that 
accreted onto GALEX1931 to be ^ "[^^^^"^ <0.2 or <1 % of the mass of the parent body. It is 



noted that use of the lower Mg abundance derived by IVennes et al\ (120111 ) would result in a 
more water-rich body. Thus it appears as though the material accreting onto GALEX1931 
is very likely to have originated in rocky minerals and is likely to be only a small percentage 
of water by mass. 

Examination of the bulk elemental composition of the material that is accreted by 
GALEX1931 can potentially yield insight into the evolutionary processes that affect terrestrial- 
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like planets during the late stages of stellar evolution. In Figure [6] we plot the elemental 
abundances rel ative to Fe for thr ee white dwarfs: GD 362 ( iZuckerman et a/.l 120071 : iKoester 
2009h . GD 40 falein et a/.ll201oh . and GALEX1931. We normalize each of these relative 
abundance measurements (at the time of accretion for each object ) to elemental abundances 
relative to Fe for CI Chondrites (see Table [3] and lloddersl l2003h . CI Chondrites are cho- 
sen as these objects have some of the best-determined bulk compositions of any object in 
our Solar system. CI Chondrites have near-Solar composition s and thus rep resent primitive 
rocks that likely formed directly from the Solar nebula (e.g., iLodderd l2003l and references 
therein). Under the assumption that primitive rocks are similar in all planetary systems, we 
seek to probe whether or not the parent body being accreted by GALEX1931 has undergone 
significant processing by comparing the GALEX1931 measured abundances to CI Chondrite 
abundances. To further place the white dwarf measurements into context, w e overplot the 



bulk e lemental composition o f Earth as extracted fron i the measurements of lAUegre et al. 
(119951 ) and the compilation of Ivan Thienen et al\ (120071 ). 



Most striking in Figure [6] is the rough agreement of the abundances of all elements 
for the three white dwarfs with the bulk Earth, especially the depletion of the volatiles C 
and O relative to what would be expected from the more primitive CI Chondrites. Next 
is the good agreement between transition metal abundances (Mn, Cr, Fe, and Ni) for the 
accreted objects and the bulk Earth. The obvious discrepancies are with Si and Mg and the 
refractories Ca and Ti. Assuming the parent bodies that are accreted by each white dwarf 
started close to Earth-like in bulk elemental composition (a risky assumption, but the best 
that can be made at present), significant deficits of Si and varying degrees of deficient Mg 
suggest crustal and mantle loss. This result suggests that each white dwarf has accreted 
objects massive enough to experience differentiation. That is, the atoms layered themselves 
in the parent body by mass (with the heaviest atoms at the parent body's core). Such a 
result indicates that the parent bodies now being accreted by each white dwarf were indeed 
massive bodies and perhaps even full-fledged planets. 

Ca, Ti, and Al abundar ices could potentia lly bear the signature of a particular crust 
and mantle loss mechanism. iKlein et al\ (120101) compar ed the accreted abundances for GD 
40 to a model developed by iFegley fc CameronI (119871 ) that described the evolution of a 
rocky planetary body when it is exposed to intense radiation. The predictions of this model 
are in good agree ment with the parent body abundances for the object accreted by GD 40 
(IKlein et a/.ll2010l ). Such a model could potentially provide a reasonable explanation for the 
parent body composition of the object accre ted by GD 362 as it too sh ows the characteristic 
enhancement of Ca and Ti predicted in the iFegley fc Cameron! (119871 ) model, although the 
mantle loss appears to be more severe as per the significant Mg deficiency (Figure [6]). To 
explain the results for GALEX1931 requires either additional evolutionary processes or a 
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different evolutionary path that results in crust and mantle loss without the enhancement 
of Ca and Ti. Indeed, GALEX1931 shows a significant deficit of Ca, Ti, and Al relative to 
the Earth and CI Chondrites (see Figure E]). 

To explain the distinctive elemental mixture accreting onto GALEX1931 we propose 
a scenario in which the dominant evolutionary mechanism acting on the parent body was 
interaction with the host star stellar wind. After irradiation by the luminous host star while 
on the asymptoptic giant branch (AGB), interaction with the A GB ejecta seems a likely 



choice for a rocky planet mass loss mechanism (e.g., I Jural 120081 ). In such a scenario the 
rocky planet is essentially sanded down by the AGB wind, losing first its crust then mantle 
depending on the severity of the planet's interaction with the stellar wind. We estimate to 
first order the relative abundances of such an object after crustal and mantle stripping. We 
begin with the assumption that the object started with the s ame composition as the bulk 



Earth and the same crust and mantle abundances as the Earth (j Anderson! 1 1 9 89l : lAUegre et al. 



19951 ). Also inherent is the assumption that the parent body's tidally shredded remnants are 
well mixed in the disk. We remove varying fractions of the crust, then upper mantle, and then 
lower mantle. After this removal step, the bulk abundances of the object are recomputed 
relative to Fe. Experimentation with varying degrees of crust and mantle removal suggests 
that the best match between the simple stripping model and the GALEX1931 abundances 
shown in Figure [6] is obtained after stripping the entire crust and upper mantle and half of 
the lower mantle. The simple model provides good agreement with the measurements, with 
the exception of aluminum. These results suggest that "wind stripping" of the outer layers 
of an Earth-like planet is a viable mode to produce the observed parent body abundances 
in at least some white dwarfs and thus that interaction with the AGB stellar wind is likely 
an important process in a terrestrial planet's final evolutionary phase. Exploration of the 
relevant physics in the wind stripping process can further develop the simple model laid out 
above. 



5. Conclusions 

We have carried out a comprehensive suite of observations aimed at characterizing 
GALEX1931's infrared excess and atmospheric pollution. Our broad data set allows us 
to identify the excess infrared emission detected toward GALEX1931 as emanating from an 
orbiting dusty debris disk. Analysis of the heavy metals detected in GALEX1931's atmo- 
sphere suggests that it is accreting the remnants of a differentiated — and hence massive 
— rocky body. Assuming the object being accreted by GALEX1931 started out roughly 
Earth-like in composition, comparison of the distinctive bulk elemental composition of this 
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terrestrial-like planet to predictions from a simple wind-stripping model suggests that the 
planet's outer layers were removed before the planet was accreted by the white dwarf star. 
A likely culprit for such stripping is interaction between the planet and its host star's dense 
stellar wind while the star progresses through the asymptoptic giant branch. 
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Fig. 1. — Images of GALEX1931 (center star in each frame). North is up and East is 
left in each image. Top row: BVRI (from left to right) images from the Nickel 1-m. All 
images are presented with the same inverted linear stretch. Bottom row: JgHRgL' (from left 
to right) images from ISAAC. The JsHKg images are all presented with the same inverted 
linear stretch. Note that the two objects within ?a2" of GALEX1931 are red. 



- 18 - 




Fig. 2. — Best fit models (red lines) to Balmer lines in the KAST blue side data of 
GALEX1931 (underlying black spectrum). The best fit Tcs and log^f are also displayed. 



1 



2 



3 



4 




Fig. 3. — Spectral energy distribution for the DAZ white dwarf GALEX1931. Overplotted 
are models for an atmosphere of a Teff=23,500 K, \og{g)=8.0, DAZ white dwarf with an 
orbiting fiat, passive, opaque dust disk having Tinnor=1400 K, Touter=1200 K, and i=70°. 
Data points are GALEX NUV and FUV (Section 13. II discusses the disc repancy between the 
model and the GALEX NUV and FUV measurements). Nickel BRI, the lVennes et all koid ) 
V, and ISAAC JgHKsL'. The red curve in the optical is the KAST spectrum. The red 
curve in the infrared is the FIRE prism-mode spectrum. Both the KAST and FIRE data 
are flux calibrated to the broadband photometry. The blue, dotted curve is the white dwarf 
atmospheric model and the blue dashed curve is the disk model. The sum of the two models 
is shown with the green dash-dotted line. 
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Fig. 4 . — Elements identified in the HIRES spectra of GALEX1931 not seen by lVennes et al. 



teOlOl ). Aluminum is tentatively identified, and the indicated abundance is suggested as an 
upper limit. Higher S/N spectra in this region or ultraviolet spectra could confirm the 
detection of Al. Wavelengths are in the vacuum heliocentric reference frame. 
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Fig. 5. — Comparison of near-infrared measurements to models for a roughly mid-L type 
brown dwarf orbiting GALEX1931. The red, upper data points are the 2MASS photome- 
try; the black, lower data points are the ISAAC photometry. The black dotted curves at 
the bot tom of the plot are mid-L type brown dwarf spectral data from the SpeX spectral 
library (ICushing et a/.ll2005l : iRayner et a/.l 120091 ) shifted in flux to the distance estimated for 
GALEX1931 (Section 13. ip . The dotted curve in the middle of the plot is an atmospheric 
model for GALEX1931. The blue and gold curves are the sum of the white dwarf and brown 
dwarf curves for an L5 and L4.5 type brown dwarf, respectively. The FIRE data smoothly 
connect the ISAAC data points. 
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Fig. 6. — Comparison of wliite d warf abundances r elative to Fe to CI Cliondrite abundances 
relative to Fe (after Figure 15 of iKlein et a/.ll2010l ). Elements have increasing condensation 
temperature to the right. The triangle data points are measurements and upper limits for 
GALE X1931 from Table [3l The diamond data points are measurements and upper limits for 

le circle data points are measurements and upper limits for 
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20071: 



GD 362 JZu ckerman 
jAUegre et a/..,1995i : 
model as described in Section HJ 
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van Thienen et al. 



20091) . The dashed line is data for the bulk Earth 



20071 ). The dotted line is a "wind-stripped Earth'' 
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Table 1. Broad-band Fluxes for GALEX1931 



Band A Fobs 
(nm) (mag) (mJy) 



ISAAC- NIR' 



L' 


3780 


14.00±0.14 


0.63±0.08^ 


Ks 


2160 


14.62±0.05 


0.95±0.05 


H 


1650 


14.80±0.05 


1.26±0.06 


Js 


1240 


14.80±0.05 


1.96±0.10 


Nickel 40'"^ - Optical 


I 


798.2 


14.47±0.04 


3.94±0.16 


R 


640.7 


14.31±0.03 


5.79±0.17 


V 


544.8 


14.20±0.02 


7.52±0.15 


B 


436.3 


14.18±0.02 


8.55±0.17 


GALEX^ - Ultraviolet 


NUV 


227.1 


13.70±0.03 


11.94±0.36 


FUV 


152.8 


13.46±0.05 


14.99±0.75 



'^ISAAC and Nickel magnitudes are on 
the Vega system. GALEX measurements 
are in AB magnitu des and have uncertain - 
ties as suggested in lMorrissey et a/.l (120071 ). 

''This value disagrees with the WISE 
3.35 Atm measured value of 1.02±0.04 mJy 
(jPebes et a/.ll2010l ). see Section [3^ 



Table 2. Spectroscopic Observations Summary 



UT Date 


Instrument 


Setup 


Coverage 


Resolution*^ 


Integration Time (s) 


S/N 


A of S/N^ (A) 


05 July 2010 


HIRES 


UV Collimator 


3130-5940 A 


~40,000 


2x1800 


80^= 


3450 


19 September 2010 


FIRE 


Prism-mode 


0.7-2.5 /im 


^250-350^^ 


4x60 


80^^ 


10500 


09 October 2010^ 


KAST 


600/4310 
1200/5000 


3430-5510 A 
5520-6910 A 


~3.4 A 
~2.4 A 


600 
600 


100 
60 


4500 
6300 



'^Resolutions were measured from the FWHM of single arclines in our comparison spectra. 
Wavelength where S/N measurement was made in the spectrum. 
'^S/N for combined exposures. 



to 



'^See iBurgasser et al\ (120101 ). 

^Observations used the 2.0" slit and the d55 dichroic. 



Table 3. GALEX1931+0117 Metal Pollution 



z 


l0g[ri(Z)/n(H)Jincasured 


Tdiff (days)'" 


[n(Z)/ra(Fe)Jaccreted 


[n(Z)/n(Fe)]ci" 


Macc/(10^ g 






o.oi 1 


(171 9 


n ^Qi 9 

U.Oc/ 


^1 1 ^ 




o.uoznu. xu 




-L . 1 OOO 


O. 1 Uc/U 


1 09 n 


A/To- 
iVig 


— zL 1 n+n 1 n 




n ^1/in 


1 9099 


4fi 1 


Al 


<-5.85 


5.336 


<0.0077 


0.0954 


<0.772 


Si 


-4.35±0.11 


4.542 


0.2884 


1.1481 


29.8 


Ca 


-5.83±0.10 


2.730 


0.0158 


0.0691 


2.34 


Ti 


<-7.00 


2.960 


<0.0010 


0.0027 


<0.174 


Cr 


-5.92±0.14 


2.560 


0.0137 


0.0151 


2.55 


Mn 


-6.26±0.15 


2.292 


0.0070 


0.0104 


2.64 


Fe 


-4.10±0.10 


2.329 


1.00 


1.00 


205.0 


Ni 


<-5.60 


1.982 


<0.0371 


0.0549 


<8.04 



^"Diffusion constants, see iKoested (120091 ) and Section HI 

'^Parent body abundances relative to Fe; see Section [Hand Figure [61 



'^CI chondrite data from iLodderj (120031 ). See also Figure [61 

'^Mi,cc{Z)=Menv{Z)/Tdis{Z) where M^nviZ) is the mass of eleme nt Z in GALE X1931's envelope 
assuming the hydrogen-dominated envelope mass is 9.4 x 10^^ g ( jKoesterl [20091 D. Koester 2010 
private communication) . 



^This value is derived from the equivalent widths reported in [Vennes et al\ ( 12010[ ). see Section [XH 



